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ABSTRACT 

We present metallicity gradients in 49 local field star-forming galaxies. We derive gas- 
phase oxygen abundances using two widely adopted metallicity calibrations based on the 
[Oni]/H^, [N Il]/Ha and [Nll]/[Oll] line ratios. The two derived metallicity gradients are 
usually in good agreement within ±0.14 dex (R 25 is the B-band iso-photoal radius), 
but the metallicity gradients can differ significantly when the ionisation parameters change 
systematically with radius. We investigate the metallicity gradients as a function of stellar 
mass (8 < \og{lVR/ Mq) < 11) and absolute B-band luminosity (—16 > Mb > —22). 
When the metallicity gradients are expressed in dex kpc“^, we show that galaxies with lower 
mass and luminosity, on average, have steeper metallicity gradients. When the metallicity 
gradients are expressed in dex R 25 , we find no correlation between the metallicity gradients, 
and stellar mass and luminosity. We provide a local benchmark metallicity gradient of field 
star-forming galaxies useful for comparison with studies at high redshifts. We investigate the 
origin of the local benchmark gradient using simple chemical evolution models and observed 
gas and stellar surface density profiles in nearby field spiral galaxies. Our models suggest 
that the local benchmark gradient is a direct result of the coevolution of gas and stellar disk 
under virtually closed-box chemical evolution when the stellar-to-gas mass ratio becomes 
high (2^ 0.3). These models imply low current mass accretion rates (< 0.3 x SFR), and low 
mass outflow rates (< 3 x SFR) in local field star-forming galaxies. 
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1 INTRODUCTION 

The content of heavy elements in a galaxy is one of the key 
properties for understanding its formation and evolutionary his¬ 
tory. The gas-phase oxygen abundance in the interstellar medium 
(ISM) of a galaxy (or “metallicities”), defined as the number 
ratio of oxygen to hydrogen atom and commonly expressed as 
12 + log(0/H), is regulated by various processes during the evo¬ 
lutionary history of a galaxy. While the oxygen is predominately 
synthesised in high-mass stars (> 8Mq) and subsequently released 
to the ISM by stellar winds and supernova explosion, the oxygen 
in the ISM could also be expelled to the circumgalactic medium, 
or potentially become gravitationally unbound, via feedback pro¬ 
cesses, e.g. galactic-scale outflows. Gas inflows triggered by merg¬ 
ers and inflows of pristine gas from the intergalactic medium could 
also dilute the metallicity of a galaxy (e.g. Kewley et al. 2006a; 


Rupke et al. 2008; Kewley et al. 2010; Rupke et al. 2010a,b). The 
relations between metallicity and other fundamental properties of 
galaxies can place tight constraints on the processes governing the 
evolution of galaxies. 

The correlation between global metallicity and stellar mass 
in star-forming galaxies, i.e. the mass-metallicity relation, is one 
of the fundamental relations for measuring the chemical evolu¬ 
tion of galaxies (Lequeux et al. 1979; Tremonti et al. 2004). Whilst 
the mass-metallicity relation was first established locally, mod¬ 
ern spectroscopic surveys have enabled precise measurements of 
metallicity out to high redshifts on large numbers of galaxies 
(e.g., Savaglio et al. 2005; Erb et al. 2006; Zahid et al. 2011, 2013, 
2014b; Wuyts et al. 2014; Steidel et al. 2014; Sanders et al. 2014) 
and laid the foundation for subsequent investigation into the phys¬ 
ical origin of the relation. Various physical processes including 
metal enriched outflows (e.g., Larson 1974; Tremonti et al. 2004), 
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accretion of metal-free gas (e.g. Dalcanton et al. 2004), and varia¬ 
tion in the initial mass function (Koppen et al. 2007) or star forma¬ 
tion efficiency (e.g., Brooks et al. 2007; Calura et al. 2009) have all 
been proposed to be responsible for shaping the mass-metallicity 
relation. The mass-metallicity relation may also have an additional 
dependency on star formation rate (SFR; e.g., Mannucci et al. 
2010; Lara-Lopez et al. 2010; Yates et al. 2012). Spatially resolved 
studies have shown that the mass-metallicity relation also holds 
on smaller scales for individual star-forming regions within galax¬ 
ies (Rosales-Ortega et al. 2012). Recent work suggests that the 
mass-metallicity relation could be a direct result of some more 
fundamental relations between metallicity, stellar and gas content 
(Zahid et al. 2014c; Ascasibar et al. 2014). 

The spatial distribution of metals in a disk galaxy can also 
provide critical insight into its mass assembly history. Disk galax¬ 
ies in the local Universe universally exhibit negative metallic¬ 
ity gradients, i.e., the centre of a galaxy has a higher metallicity 
than the outskirts (e.g., Zaritsky et al. 1994; Moustakas et al. 2010; 
Rupke et al. 2010b; Sanchez et al. 2014, and references therein). In 
several cases, where measurements are possible out to very larger 
radii (2 x R 25 ^), the metallicity gradients flatten in the outer disks, 
suggesting inner-to-outer transportation of metals via mechanisms 
such as galactic fountains (e.g., Werketal. 2011; Bresolin et al. 
2012; Kudritzki et al. 2014; Sanchez et al. 2014). Extreme exam¬ 
ples of metal mixing occurs in interacting galaxies, where the non 
axis-symmetric potential induces radial inflows of gas. Both ob¬ 
servations and simulations confirm that mergers of disk galaxies 
present shallower metallicity gradients than isolated disk galaxies 
due to effective gas mixing (e.g., Kewley et al. 2010; Rupke et al. 
2008, 2010a,b; Torrey et al. 2012; Rich et al. 2012). 

Sophisticated modelling of the evolution of metallicity gra¬ 
dients in disk galaxies has shed light on the formation, gas ac¬ 
cretion, and star formation history of the disks. While the details 
vary from model to model, typical assumptions of inside-out disk 
growth, no radial matter exchange, and closed-box chemical evolu¬ 
tion can successfully reproduce the current gradients in local galax¬ 
ies (e.g., Chiappini et al. 2001; Fuetal. 2009). However, some 
models predict that metallicity gradients steepen with time (e.g., 
Chiappini et al. 1997, 2001; see also Mott et al. 2013 who included 
radial inflow), while others predict the opposite (e.g., Molla et al. 
1997; Prantzos & Boissier 2000; Fuetal. 2009; Pilkington et al. 
2012). Testing the model predictions using observations of high 
redshift galaxies are challenging (e.g.. Yuan et al. 2011; Jones et al. 
2010, 2013, and references therein). In addition, systematic effects 
from insufficient resolution and/or binning of the data unfortunately 
can seriously the reliability of metallicity gradients measured at 
high redshifts (Yuan et al. 2013; Mast et al. 2014). 

Statistical studies of metallicity gradients in the local Uni¬ 
verse provide an alternative approach to constrain the theoreti¬ 
cal simulations. Although the measurements are time-consuming, 
sample sizes of few tens of galaxies have been achieved in the 
past using long-slit spectroscopy. These samples gave intriguing 
(and sometimes contradictory) correlations between metallicity 
gradients and physical properties of the disk galaxies. For exam¬ 
ple, barred galaxies tend to exhibit shallower metallicity gradi¬ 
ents than non-barred galaxies even when galaxy sizes are taken 
into account (e.g., Vila-Costas & Edmunds 1992; Zaritsky et al. 
1994), but such discrepancy is insignificant in some recent studies 
(Sanchez et al. 2014). For unbarred galaxies, galaxies with higher 

^ Radius of the 25th magnitude/arcsec^ isophote in B-band. 


6-band luminosity or higher total mass have shallower metallic¬ 
ity gradients (Vila-Costas & Edmunds 1992; Garnett etal. 1997); 
however, such behaviour is not pronounced in some measure¬ 
ments (e.g., van Zee et al. 1998; Prantzos & Boissier 2000). Some 
studies find that non-barred galaxies show a statistically signif¬ 
icant correlation between metallicity gradient and Hubble Type, 
where earlier types have shallower metallicity gradients than later 
types (Vila-Costas & Edmunds 1992; Oey & Kennicutt 1993), but 
considerable scatter exists in other measurements (Zaritsky et al. 
1994). Most studies find no correlation when metallicity gradients 
are normalised to some scale-length (i.e., R 25 , the disk scale-length 
Rd, or the effective radius Re^)- The contradictory results of some 
earlier studies might be due to the small sample sizes and incon¬ 
sistent methodologies of measuring metallicity gradients. Applying 
different metallicity diagnostics can introduce considerable system¬ 
atic errors (Kewley & Ellison 2008). 

Advances in instrumentation such as multi-slit spectroscopy 
and wide-field integral field spectroscopy (IFS) is in the pro¬ 
cess of revolutionising statistical studies of metallicity gradients 
in the local Universe (e.g., Sanchez et al. 2014). Large on-going 
and future large IFS surveys include the Calar Alto Legacy In¬ 
tegral Field Area Survey (CALIFA; Sanchez et al. 2012a), the 
Sydney-AAO Multi-object Integral field spectrograph (SAMI) Sur¬ 
vey (Groom et al. 2012; Bryant et al. 2014; Allen et al. 2014), the 
Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) 
Survey, the Hector Survey (Lawrence et al. 2012; Bland-Hawthom 
2014), and many others. These IFS surveys are not only extremely 
efficient in collecting large numbers of spectra simultaneously and 
seamlessly across an entire galaxy, but also have desirable wave¬ 
length coverage to capture multiple key emission lines for deriving 
metallicity. Such features pose a unique opportunity to eliminate 
systematic errors using statistical approaches. 

In this paper, we study the metallicity gradients in a sample of 
49 local field star-forming galaxies. We derive metallicity gradients 
using different abundance calibrations and discuss potential sys¬ 
tematic effects induced by the calibrations adopted. We investigate 
whether metallicity gradients in field star-forming galaxies corre¬ 
late with their physical properties. We show that there is a common 
metallicity gradient in local field star-forming galaxies and we pro¬ 
vide some benchmark values. Finally, we adopt simple chemical 
evolution models to explain the formation of the common metallic¬ 
ity gradient. 

The paper is structured as follows. We describe our samples, 
observations and data reduction in Section 2, and our data analysis 
in Section 3. In Section 4, we detail our methodology of deriving 
the metallicity, ionisation parameter, and metallicity gradients. In 
Section 5, we present our measurements of metallicity gradients, 
discuss the systematic effects, and compare the metallicity gradi¬ 
ents with stellar masses and absolute 6-band magnitudes. We pro¬ 
vide a benchmark metallicity gradient in Section 6 and investigate 
the origin of the benchmark gradient in Section 7 using the simple 
chemical evolution models. Finally, a summary and conclusions are 
given in Section 8. Through out this paper, we assume the standard 
A cold dark matter cosmology with Hq — 70 km Mpc“^, 
Dm = 0.3 and Ha = 0.7. 


^ The effective radius is the radius at which the integrated flux is half of the 
total one. Comparing to the disk scale-length for the classical exponential 

profile. Re = 1.678356d- 
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2 SAMPLES 

The 49 galaxies studied in this work are drawn from various 
sources, including literature data, public data and our targeted ob¬ 
servations. We describe the four samples in the next four subsec¬ 
tions. The focus of this paper is to investigate field star-forming 
galaxies, and therefore we select only field galaxies that are not un¬ 
dergoing major mergers and not in close pairs; none of our galaxies 
have massive companions (i.e. stellar mass higher than one-third of 
the main galaxies) within 70 kpc in projection and 1000 km in 
line-of-sight velocity separation. 


2.1 CALIFA galaxies in Data Release 1 

We use the publicly available IPS data from the first data re¬ 
lease (DRl) of the CALIFA survey. A full description of the sur¬ 
vey design, including details of target selection and data reduction 
scheme, can be found in Sanchez et al. (2012a). Specific details re¬ 
lated to the DRl can be found in Husemann et al. (2013, see also 
Walcher et al. 2014). Below we briefly summarise the information 
relevant to this work. 

The CALIFA DRl contains reduced IPS data of 100 local 
galaxies (0.005 < 2 < 0.03) covered by the Sloan Digital Sky Sur¬ 
vey (SDSS; York et al. 2000; Abazajian et al. 2009). In this study, 
we focus only on spiral galaxies that are not undergoing a major 
merger, and not in close pairs. Extremely edge-on systems with in¬ 
clination angle larger than 70° are excluded from our analysis since 
de-projecting radial distance is uncertain. Systems without enough 
sufficiently high signal-to-noise spaxels (S/N > 3) to measure 
emission line fluxes are also excluded. In total, 21 CALIFA galax¬ 
ies are analysed. 

The released CALIFA datacubes, as processed by the CAL¬ 
IFA automatic data reduction pipeline, have a spatial size of ~ 
74" X 64" on a rectangular 1" grid. The point spread func¬ 
tion, as measured from field stars in the datacubes, has a me¬ 
dian full-width measured at half-maximum (FWHM) of 3.7". Ev¬ 
ery CALIEA galaxy is observed using the fibre bundle integral 
field unit (IFU) PPak, and two different setups with the Pots¬ 
dam Multi-Aperture Spectrophotometer (PMAS; Roth et al. 2005; 
Kelz et al. 2006). The low-resolution (V500) and high-resolution 
(V1200) setup cover wavelength ranges of ~ 3745-7500 A and 
~ 3650-4840 A, respectively. The V500 reduced datacubes have a 
FWHM spectral resolution of 6.0A (R ~ 850) and a spectral chan¬ 
nel width of 2.0A. The V1200 reduced datacubes have a FWHM 
spectral resolution of 2.3A (R ~ 1650) and a spectral channel 
width of O. 7 A. 


2.2 WiFeS galaxies 

The CALIFA sample was selected based on the angular iso-photal 
diameter of the galaxies (45" < D 25 < 80"). Therefore, the 
CALIFA sample is inevitably biased towards galaxies of higher 
mass (> 10® Mq). To probe metallicity gradients in galaxies of 
lower mass in a statistically significant way, we conducted supple¬ 
mental IFS observations to specifically target low mass systems 
(i.e., log{Mt,/ Mq) = 8-9). We first selected a mother sample 
of low-mass galaxies from the SDSS Data Release 7 value-added 
catalog constructed by the MPA/JHU group®. We used the stellar 


® http://www.mpa-garching.mpg.de/SDSS/DR7/ 


mass derived by the MPA/JHU group as a reference for target se¬ 
lection (Kauffmann et al. 2003a; Salim et al. 2007); the final stel¬ 
lar masses adopted in this work are derived separately and con¬ 
sistently for all our samples. As a result of this selection, one of 
the WiFeS galaxies presented has a substantially larger final stel¬ 
lar mass (log(M*/MQ) ~ 10.2) due to the incorrect apertures 
adopted by MPA/JHU. We remove galaxies with AGN from our 
mother sample using the optical line ratios [N Il]/Ha and [O III]/H/J 
(Kewley et al. 2006b). We further constrained the mother sample 
to have low inclination disks and spatial extent comparable to the 
Wide Field Spectrograph (see Section 2.2.1). We also visually con¬ 
firmed that these galaxies are not undergoing major merger and do 
not have massive companions. From the mother sample, we then 
selected our final targets based on observability and instrumental 
sensitivity. In total, we observed 19 galaxies, 10 of which yield re¬ 
liable metallicity gradients and are presented in this paper. 



WiFeS central 3" 03N2 index 


Figure 1. Comparison of the observed 03N2 index (Equation 1) between 
those from the WiFeS and the SDSS data. Each dot corresponds to one of 
the 10 WiFeS galaxies that were also observed by the SDSS spectroscopic 
survey. To derive the 03N2 index from the WiFeS data, we extract line 
fluxes in 3"apertures at the locations of the SDSS fibres. The 03N2 index 
derived from the two datasets are consistent within approximately 0.1 dex. 


2.2.1 Obsenation and data reduction 

We observed our low-mass galaxies using the WiFeS on the 2.3- 
m telescope at Siding Spring Observatory in December 2012 and 
April 2013. WiFeS is a dual beam, image-slicing IFU consisting 
of 25 slitlets. Each slitlet is 38" long and 1" wide, yielding a 
25" X 38" field of view. For a thorough description of the instru¬ 
ment, see Dopita et al. (2007) and Dopita et al. (2010). Each galaxy 
was observed with the blue and red arms simultaneously using 
the B3000 and R7000 gratings, respectively. All galaxies were ob¬ 
served with a single WiFeS pointing except for 1031752.75-071804 
where we adopted a two-point mosaic. The typical exposure time 
is ~ 1 — 2 hours per galaxy under seeing conditions of 1.5 — 2.5". 

We reduce the data using the custom-built data reduction 
pipeline PYWIFES (Childress et al. 2014). The final reduced data 
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consist of two datacubes on 1" x 1" spatial grids for each galaxy. 
The blue cube covers ~ 3500-5700A with a FWHM velocity reso¬ 
lution of ~ 100 km at H/3 1.7A or R ~ 3000) and a spec¬ 

tral channel width of ~ O.SA. The red cube covers ~ 5500-7000A 
with a FWHM velocity resolution of ~ 40 km at Her (~ O.OA 
or R ~ 7000) and a spectral channel width of ~0.4A. 

To compare our reduced WiFeS datacubes with the SDSS fibre 
spectra, we present in Figure 1 the emission line ratios, the 03N2 
indexes (see below; Equation 1), derived from the two datasets. For 
the SDSS data, we adopt the line fluxes from the MPA/JHU value- 
added catalog; for the WiFeS data, we extract line fluxes within the 
corresponding fibre apertures from the emission line maps (see be¬ 
low). Figure 1 demonstrates that the 03N2 indexes from the WiFeS 
and SDSS data are consistent within approximately 0.1 dex. 

2.3 Galaxies from Sanchez et al. (2012b) 

To further increase our sample size, we analyse 9 field galaxies pre¬ 
sented in Sanchez et al. (2012b, hereafter S12) with low inclination 
angles « 50°). 

S12 studied ~ 2600 Hit regions in 38 nearby galaxies se¬ 
lected from the PINGS survey (Rosales-Ortega et al. 2010) and 
Marmol-Queralto et al. (2011). All 38 galaxies were observed with 
PPak and PM AS, which deliver spectral coverages similar to the 
CALIFA data of ~ 3700 — 6900A. S12 applied a semi-automatic 
procedure, HiiExplorer, to search for Hit regions in lEU data 
under the assumptions that Hit regions are peaky and isolated 
line-emitting structures with typical physical sizes of few hun¬ 
dred parsecs (more details in S12). A very similar spectral fit¬ 
ting approach (to Section 3.1) was applied on synthetic spectra of 
Hit regions to decouple the underlying stellar contribution from 
line emissions. The final public flux catalogs contain seven strong 
lines, [On] AA3726,3729, H^S, [Out] A5007, [Oi] A6300, Ha, 
[Nil] A6583, and [Sii] AA6716,6731. 

All the 9 galaxies analysed in this study have multiple 
bright Hll regions measured in all the strong lines including 
[O II] AA3726,3729, which allows reliable constraints simultane¬ 
ously on metallicity and ionisation parameters with different diag¬ 
nostics. 

2.4 Galaxies from Rupke et al. (2010b) 

Rupke et al. (2010b, hereafter RIO) measured metallicity gradients 
in interacting and non-interacting galaxies. They show that, on av¬ 
erage, interacting systems present shallower metallicity gradients 
than non-interacting systems. Their control sample comprises 11 
non-interacting local galaxies and they measure their metallicity 
gradients using published emission line data from Hll regions. Two 
of these control sample galaxies overlaps with the S12 sample. 
We include their measurements for the remaining 9 galaxies in 
our analysis. Their methods of correcting for extinction and de¬ 
riving metallicity are exactly the same as our work and therefore 
we include their measurements of metallicity gradients without a 
re-analysis of the data. 


3 ANALYSIS 

3.1 Emission line maps 

To place constraints on metallicity and ionisation parameter using 
emission line diagnostics, we measure emission line fluxes in each 


spaxel of CALIEA and WiEeS galaxies by spectral fitting. We use 
an earlier version of the spectral fitting tool LZIFU described in 
Ho et al. (in preparation; see also Ho et al. 2014). The fitting ap¬ 
proaches for the two samples are very similar, though some details 
are fine tuned to accommodate the differences in spectral coverage 
and resolution between the two datasets. Below, we first elaborate 
our method for the CALIEA galaxies before describing the different 
treatments for the WiEeS galaxies. 

Prior to fitting the CALIFA galaxies, we first correct for the 
known spatial misalignment between the V500 and VI200 dat¬ 
acubes (Husemann et al. 2013). For each CALIFA galaxy, we re¬ 
align the two cubes by correlating continuum images constructed 
separately from the two cubes using a wavelength range covered 
by both settings (4240-4620A). Typical misalignments are ~ 1” 
to 2" while in several extreme cases ~ 3” to 5”. We then rescale 
the VI200 data to match the V500 data using scale factors deter¬ 
mined from the median flux density ratios between an overlapped 
spectral coverage of 4000A and 4500A. Slight differences in flux 
levels may be a result of imperfect calibration and any residual sys¬ 
tematic errors. 

We then perform simple stellar population (SSP) synthesis to 
remove the underlying stellar continuum before fitting emission 
lines. To model the stellar continuum, we adopt the penalised pixel¬ 
fitting routine (PPXF; Cappellari & Emsellem 2004) and employ 
the empirical MIUSCAT'* SSP models of 13 ages® and 4 metallic- 
ities® (Vazdekis et al. 2010) while assuming a Salpeter initial mass 
function (IMF; Salpeter 1955). Bad pixels, sky lines and the vicin¬ 
ity (ilSA) of emission lines are masked prior to fitting the con¬ 
tinuum models. After subtracting the continuum models from the 
data, we further remove low frequency fluctuations in the residuals 
by fitting fourth order b-spline models. We note that the major goal 
of performing SSP synthesis fitting is to correct for stellar Balmer 
absorption; we do not derive stellar age and metallicity from our 
SSP fits. 

Subsequent to removing the continuum in the V500 
and VI200 datacubes separately, we model the strong emis¬ 
sion lines ([On] AA3726,3729, H,5, [Oiii] AA4959,5007, 
[Nil] AA6548,6583, Ha, and [S II] AA6716,6731) as simple gaus- 
sians. We perform a bounded value nonlinear least-squares fit using 
the Levenberg-Marquardt method implemented in IDL (Markwardt 
2009, MPFIT). We constrain (1) all lines to have the same velocity 
and velocity dispersion, (2) the ratios [N II] A6583/[N II] A6548 
and [O III] A5007/[0 III] A4959 to their theoretical values given by 
quantum mechanics, (3) the velocity to be between -1-600 km 
and —600 km to the systemic velocities as measured from 
SDSS, and (4) the velocity dispersion to be between 50 and 1000 
km s“^. 

Figure 2 shows a typical spectrum and spectral fit. This fig¬ 
ure best illustrates the above procedure of decoupling stellar con¬ 
tribution in each spaxel from emission lines originated predomi¬ 
nantly from Hll regions. Figure 3 shows two emission line maps, 
SDSS composite image, extinction map, velocity field map, and a 
key diagnostics line-ratio map of the CALIFA galaxy NGC7321 to 
demonstrate the final products from our analysis described above. 

For fitting the WiFeS galaxies, the above procedure is adopted 
with some minor modifications to accommodate the different spec- 


^ http://miles.iac.es/pages/ssp-models.php 

® 0.06, 0.10, 0.16, 0.25, 0.40, 0.63, 1.00, 1.58, 2.51, 3.98, 6.31, 10.00, and 
15.85 Gyr 

® [M/H] = -0.71,-0.40, 0.00, and 0.22 
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Figure 2. An example of the spectral fitting approach applied on the CALIFA data (see Section 3.1 for details). The grey and black thick lines in the upper 
panel indicate the wavelength ranges where the data (black: V500 data; grey: V1200 data) are adopted to constrain the continuum models (red: V500; pink: 
V1200). Bad channels, the vicinity of strong emission lines and sky lines are excluded from the fit. The four middle panels show emission lines (red) fit to the 
continuum subtracted spectra (black). All lines are fit simultaneously and share the same velocity and velocity dispersion. The bottom four panels show the 
residuals, and the blue dashed lines indicate the ztlcr noise levels. 


tral coverages and resolutions. No correction for misalignment is 
required for the WiFeS datacubes since the blue and red data were 
observed simultaneously. SSP synthesis fitting is performed simul¬ 
taneously on both cubes to take full advantage of the 4000A break, 
an important age indicator, captured only in the blue data. Fitting 
the red side separately would in principle increase the degenera¬ 
cies in SSP synthesis fitting. We first down-grade the red data to 
the same spectral resolution as the blue data (R ~ 3000), mask 
out noisy parts of the spectra due to poor CCD sensitivities, and 
merge the two datacubes to form a master datacube which covers 
~ 3700-6950A. We then use PPXF and theoretical SSP models, 
assuming Padova isochrones, of 18 ages^ and 3 metallicities® from 
Gonzalez Delgado et al. (2005) to determine the contribution from 
different stellar populations and the degree of dust extinction. The 
results are then used to reconstruct continua of the blue and red data 
at their native spectral resolution. The same line fitting approach is 
then applied to the continuum subtracted cubes and yields emission 
line maps. 

3.2 Extinction correction 

We use a consistent method for all the galaxies to correct the 
wavelength-dependent extinction caused by dust attenuation. For 

’’ 0.004, 0.006, 0.008, 0.011, 0.016, 0.022, 0.032, 0.045, 0.063, 0.089, 
0.126, 0.178, 0.251, 0.355, 0.501, 0.708, 1.000, and 1.413 Gyr. 

® Z = 0.004, 0.008, and 0.019. 


a given measurement (i.e., a spaxel for the CALIFA and WiFeS 
samples or an Hll region for the S12 samples), we assume the 
classical extinction law by Cardelli etal. (1989) with Rv = 3.1 
and Ha/H/I = 2.86 under the case-B recombination of Te = 
10, 000 K and ne = 100 cm“® (Osterbrock & Ferland 2006). This 
prescription is consistent with that adopted in RIO. 

3.3 Other physical quantities 

3.3.1 Stellar mass (M*) 

We use the Le Phare® code developed by Amouts, S. & Ilbert, 
O. to estimate the galactic stellar mass. Le Phare compares pho¬ 
tometry measurements with stellar population synthesis models, 
based on a template-fitting procedure, to determine mass-to- 
light ratios, which are then used to estimate the stellar mass of 
galaxies. The stellar templates of Bruzual & Chariot (2003) and 
a Chabrier IMF (Chabrier 2003) are used to synthesise magni¬ 
tudes. The 27 models span three metallicities and seven expo¬ 
nentially decreasing star formation models (SFR oc with 

r = 0.1, 0.3,1, 2, 3, 5,10,15 and 30 Gyr. We apply the dust atten¬ 
uation law from Calzetti et al. (2000) allowing E{B — V) to vary 
from 0 to 0.6 and stellar population ages ranging from 0 to 13 Gyr. 

Photometric measurements are collected from various 
sources. For all the 21 CALIFA galaxies, 10 WiFeS galaxies, and 
5/9 S12 galaxies, we adopt values from the SDSS DR7 photometry 

® http://www.cfht.hawaii.edu/~ amouts/LEPHARE/lephare.html 
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Figure 3. An example of 2D maps from our spectral analysis described in Section 3.1. The first row shows the Ha map, [0II] AA3726,3729 map, and SDSS 
3-colour image of NGC7321, one of the CALIFA galaxies. The second row shows the E(B-V), velocity field, and 03N2 maps. The bright foreground star in 
the SDSS image is masked out in all the other maps. 


catalog (Abazajian et al. 2009) and 2MASS extended source cata¬ 
log (Skrutskie et al. 2006). The Petrosian magnitudes of SDSS u, g, 
r, i, and z-band are corrected for the foreground Galactic-extinction 
(Schlegel et al. 1998). The 2MASS J, H, and Ks magnitudes mea¬ 
sured from fit extrapolation are adopted to approximate the total 
magnitudes. All these galaxies have the 5-band SDSS photometry 
and the majority (28/36) also have the 3-band 2MASS photometry. 
For the rest of the four S12 galaxies and the nine RIO galaxies, we 
collect available U, B, V, and the 3-band 2MASS photometric mea¬ 
surements from NASA/IPAC Extragalactic Database (NED). All 
these galaxies have the 3-band 2MASS photometry. The B and V- 
band photometry are available for all galaxies, and the D-band pho¬ 
tometry is available for about half of these galaxies (6/13). These 
optical photometric measurements are also corrected foreground 
Galactic-extinction. When reasonable uncertainties of the measure¬ 
ments are unavailable, we assume 0.1 dex for the optical bands and 
0.05 dex for the infrared bands. 

We note that the uncertainties in stellar mass are typically 
dominated by systematic errors. Different stellar mass estimators 
employ different algorithms and stellar libraries, but the estimated 
M* typically agrees within ~ 0.3 dex, marking the degree of 
systematic errors in the measurements (e.g., Drory et al. 2004; 
Conroy et al. 2009). 


3.3.2 Inclination angle 

Inclination angles of CALIFA galaxies are estimated using the con¬ 
version provided by Padilla & Strauss (2008). The effects of dust 
extinction and reddening were taken into account in their analysis. 
An estimate of inclination angle for spiral galaxy is inferred from 
the measured axis ratio {b/a) and r-band absolute magnitudes. 
Axis ratios estimated by the CALIFA team from SDSS r-band im¬ 
ages are adopted. Inclination angles of WiFeS galaxies are drawn 
from Hyperleda (Paturel et al. 2003), which assumes the classical 
Hubble formula (Hubble 1926). Inclination angles of S12 galaxies 
are taken directly from table 1 in S12, which also refers to values 
from Hyperleda. Inclination angles of RIO galaxies are taken di¬ 
rectly from table 2 in RIO, which is a compilation from various 
references. 


3.3.3 Size and distance 

We compile the sizes and distances of our samples from NED and 
Hyperleda. i ?25 from Hyperleda is adopted throughout the paper to 
quantify sizes of the galaxies. Redshift independent distances are 
available for all RIO galaxies in Table 2 of RIO. Redshift indepen¬ 
dent distances for most S12 galaxies (7/10) are adopted from NED. 
These redshift independent distances (d < 30 Mpc) are measured 
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from the Tully-Fisher relation, tip of the red-giant brach method, 
planetary nebulae, type-II supemovae, or Cepheid variables. For 
CALIFA, WiFeS, and the rest of the three S12 galaxies, we adopt 
Flubble distances inferred from their redshifts. 


4 DERIVATION OF METALLICITIES AND 
METALLICITY GRADIENTS 

In this section, we describe our methodology for using emission 
line ratios to derive metallicities, metallicity gradients and ionisa¬ 
tion parameters. 


4.1 Metallicity 

The most direct way to determine an ISM metallicity is by first 
measuring electron temperatures (Te) with temperature sensitive 
line ratios, e.g., [Olll] A4363 to [Olll] A5007, and then convert 
emission measures to metallicity after correcting for unseen stages 
of ionisation. Since [O III] A4363 is typically unavailable or only 
detected in very limited (i.e. the hottest) regions in IFU surveys, 
measuring metallicity usually relies on empirical or theoretical cal¬ 
ibrations (or a combination of both) based on strong emission lines, 
such as those available in our samples. Various such calibrations are 
available in the literatures and are widely adopted to derive metal¬ 
licity (e.g., Kewley & Dopita 2002; Pettini & Pagel 2004). 

We derive metallicities with two different calibrations elabo¬ 
rated below. Throughout the paper, we express the metallicity is in 
terms of the number ratio as 12 -|- log(0/F[). 


4.1.1 03N2 index / PP04 


The 03N2 index, defined as 


03N2 = log 


[Oiii] A5007/H/3 
[Nil] A6583/Ha ’ 


( 1 ) 


is a widely used metallicity diagnostic in the literature. An em¬ 
pirical calibration is provided by Pettini & Pagel (2004, hereafter 
PP04). The popularity of this diagnostics arises for two reasons. 
Firstly, the four lines involved are usually easily measured in lo¬ 
cal galaxies out to large radii using modern instruments. Secondly, 
because of the minimum wavelength differences between the two 
pairs of lines, the 03N2 index is virtually free from systematics 
caused by the assumption of an extinction law or reddening uncer¬ 
tainties. Nevertheless, the calibration provided by Pettini & Pagel 
(2004) is calibrated empirically with a linear fit to metallicities of 
137 extragalactic Hll regions (131 with Te-based metallicity and 6 
with detailed photoionisation models). Variation of ionisation pa¬ 
rameter, q, is not considered in the calibration. Neglect of this pa¬ 
rameter may cause serious systematic errors in the results. An up¬ 
dated calibration based on many more Te-based metallicities of Hll 
regions (309) is provided by Marino et al. (2013). This new cali¬ 
bration presents a significantly shallower slope between the 03N2 
index and metallicity than the PP04 calibration. 


4.1.2 N202 index / KD02 


The N202 index, defined as 


N202 = 


[Nii]A6583 
[On] AA3726,3729’ 


( 2 ) 


is an alternative metallicity diagnostic that has been calibrated by 
Kewley & Dopita (2002, hereafter KD02) using theoretical pho¬ 
toionisation models. We adopt the parametrisation for q = 2 x 
lO’^cm in KD02 to derive metallicity. The N202 index is in¬ 
sensitive to variation of ionisation parameter by virtue of the sim¬ 
ilar ionising potential of N+ and O^. Despite the insensitivity to 
ionisation parameter, N202 is not often used in local studies pri¬ 
marily because some spectrographs are not sensitive enough at 
~ 3700A to observe [Oil] AA3726,3729. The N202 calibration 
depends more on the assumed extinction law and extinction esti¬ 
mate than 03N2 due to the larger wavelength separation between 
[N II] A6583 and [O II] AA3726,3729. We note that the variation in 
the N/0 ratio with 0/H could be the largest uncertainty affecting 
these strong line diagnostics (Henry et al. 2000). The strength of 
the [O II] AA3726,3729 lines are strongly affected by the electron 
temperature, governed predominately by the O/H ratio and ionisa¬ 
tion parameter, while [O III] A5007 is mostly sensitive to the ioni¬ 
sation parameter. 

A long standing problem in chemical studies has been that 
different metallicity calibrations do not return consistent metal¬ 
licity measurements. Kewley & Ellison (2008) applied 10 differ¬ 
ent metallicity calibrations on the same SDSS dataset and found 
that different calibrations yield different mass-metallicity relations. 
Both the slopes and the intercepts of the mass-metallicity relation 
are significantly different from calibration to calibration. By allow¬ 
ing the mass-metallicity relations derived with different calibra¬ 
tions to be converted to the same bases, Kewley & Ellison (2008) 
derived empirical conversions between different calibrations. In 
this paper, metallicities derived using the 03N2 method are sub¬ 
sequently converted to the KD02 scale using the conversion by 
Kewley & Ellison (2008). 

We note that, when deriving metallicities using a certain diag¬ 
nostic, we only use data with S/N > 3 on all the lines associated to 
that particular diagnostic. The same rule also applies to the deriva¬ 
tion of the ionisation parameter (see Section 4.3). 

4.2 The importance of non-thermal excitation and diffuse 
ionised gas 

It is important to understand that all the above metallicity diag¬ 
nostics are calibrated empirically or theoretically using Hll re¬ 
gions. This means, by applying the diagnostics, one implicitly 
assumes that that all the nebular emission originates from pho¬ 
toionisation and heating caused by the extreme ultraviolet pho¬ 
tons emitted by O and B-type stars. If other ionisation sources 
are present, metallicity measurements would be contaminated. 
Other ionisation sources such as active galactic nuclei (AGNs) 
can have localised or even global effects on line ratios. Interstel¬ 
lar shocks originating from AGN outflows, supernovae, or stel¬ 
lar winds are also sources of non-thermal radiation and can af¬ 
fect emission line ratios. To remove measurements affected by 
non-thermal radiation from our subsequent analyses, we use line 
ratio diagnostics commonly adopted to distinguish normal from 
active galaxies (i.e., [Olll] A5007/H/3 v.s. [Nil] A6583/Ha or 
BPT diagram; Baldwin et al. 1981; Veilleux & Osterbrock 1987). 
We adopted the empirically separation line derived from SDSS 
by Kauffmann et al. (2003b, hereafter K03; see also Kewley et al. 
2001 and Kewley et al. 2006b) to exclude data contaminated by 
non-thermal excitation. 

Emission from the diffuse ionise gas (DIG; also known as 
the warm ionised medium) can be a non-negligible component 
in IFU measurements. The DIG is hot (~ 10"* K) and tenuous 
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Figure 4. An Example of determining spaxels heavily contaminated by the 
diffuse ionised gas. Details are described in Section 4.2. The black points 
are all star-forming spaxels with > Scr detections on the [SII] AA6716,6731 
and Ho lines. The blue and red squares are those spaxels with > Scr de¬ 
tections on all lines associated with the 03N2 index, i.e., [O III] A5007, 
[NII] A6583, Ho, and H/3. The green curve indicates the best fit to the 
black points using equations 3 and 4, which we adopt to determine the crit¬ 
ical Ho flux above which the covering fraction of Hll region exceeds 80%. 
Data below the critical Ho flux (red) are excluded from deriving 03N2 
metallicity. 


gas (~ 10~“ cm“®) permeating the interstellar space and extend¬ 
ing more than 1 kpc above the disk plane (see Mathis 2000 and 
Haffner et al. 2009 for reviews). The DIG generally presents much 
lower Ha surface brightness than typical Hit regions and displays 
distinctly different line ratios from Hit regions. The primary ex¬ 
citation sources are thought to be predominately ionising photons 
escaping and traveling kilo-parcsec distances from O and B-type 
stars. Measurements in an IFU spaxel can contain both emission 
from underlying Hit regions and from the DIG along the line-of- 
sight. 

To quantify the fractional contributions from the DIG and 
Hit regions in a given spaxel, we adopt a similar approach as 
Blanc et al. (2009). Spaxels dominated by the DIG are subse¬ 
quently removed from our analyses. We constrain the contribution 
from the DIG with the observed the observed [S Il]/Ha ratio 


[Sll] 

Hq 


= Z' 


^ ,[Sii], JSii], 

+ ODIG(-rr—)DIG 


Ha 


(3) 


where [S II] denotes the total flux of [S II] A6716 and [S II] A6731. 
The terms (^hh Cdig represent fractions of emission lines 
originated from Hll regions and the DIG, respectively. The sum of 
Chii and Gdig is unity. Z' denotes metallicity of the galaxy nor¬ 
malised to that of the Milky Way, i.e. Z’ = ZIZmw ■ Following 
Blanc et al. (2009), we adopt the value of ([Sll]/Ha)Hn 0.11 
and ([S Il]/Ha)DiG as 0.34. These values are supported by obser¬ 
vations of Hll regions and the DIG in our Milky Way (Madsen et al. 
2006). 

Figure 4 shows a typical [S Il]/Ha versus Ha flux plot of the 
CALIFA galaxy NGC6497. All the data points (spaxels) are sig¬ 
nificantly detected (> 3cr) in [SII] and Ha. Spaxels with high Ha 
fluxes have low (high) [Sll]/Ha, consistent with the low (high) 
line ratios of Hll regions (DIG). Hll regions are generally located 
on spiral arms and DIG is generally located in inter-arm regions. 

Similar to Blanc et al. (2009), we model the covering fractions 


of Hll regions in each galaxy as a simple function of 


Ch„ 


1 - 

/(Ha)’ 


(4) 


where /(Ha) is the Ha flux. Combining equation 3 and equation 4, 
we fit simultaneously Z' and /o, and show the best fit as the green 
curve in Figure 4. We typically find equation 4 describes the data 
reasonably well, but the theoretical basis of this functional form 
is unclear. The best fit provides a guide for imposing a criterion 
on the Ha flux to reject spaxels below a characteristic covering 
fraction of Hll regions. We exclude all spaxels below an Ha flux 
value at which the corresponding Ghu is 0.8 on the best fit curve. 
Although this criterion seems arbitrarily strict, changing the char¬ 
acteristic covering fraction to zero has only a minor impact on the 
metallicity gradients of individual galaxies, and none of our con¬ 
clusions change. 

The reason that rejecting DIG dominated spaxels has limited 
effect is that the covering fraction cut does not remove many spax¬ 
els that have not already been rejected by the S/N cuts for metal¬ 
licity diagnostics and the BPT criterion (K03) for non-thermal ex¬ 
citation. The strictest criteria in rejecting the data are the S/N cuts 
on the weak lines such as and [N II] AA6548,6583. Since the 
DIG has intrinsically low surface brightness, with the current depth 
in CALIFA and WiFeS samples, most spaxels satisfying multiple 
S/N > 3 criteria in line emissions are those with low DIG cov¬ 
ering fractions. We demonstrate this in Figure 4. The blue and red 
points are the spaxels satisfying the S/N > 3 criteria in the 03N2 
diagnostic and the criterion on the BPT (K03). The red points are 
those further rejected by the DIG criterion. 

We note that rejecting data dominated by the DIG is only re¬ 
quired for spaxel-to-spaxel analysis, i.e. the CALIFA and WiFeS 
samples. For the S12 galaxies, the fluxes were measured on ex¬ 
tracted spectra of Hll regions (i.e., binning spaxels in the vicinity 
of bright and compact Ha knots). For RIO, the flux measurements 
were performed by placing long-slits on Hll regions. DIG contam¬ 
ination is expected to be negligible in these two cases. 

In principle, after the fractional contribution is determined, 
one should be able to subtract the contribution from the DIG and 
derive the fluxes from underlying Hll regions (as in Blanc et al. 
2009). However, line ratios involving [OII] AA3726,3729 for the 
DIG are not well constrained and could produce large scatter (e.g., 
[OII] AA3726,3729/Ha; Mierkiewicz et al. 2006). Even for the 
well-measured ratios, such as [S Il]/Ha and [N Il]/Ha, considerable 
scattering and correlation between [S II]/[N II] could complicate the 
corrections and induce potential systematic errors (Madsen et al. 
2006). In this work, we adopt the simplest approach of rejecting 
inappropriate data. Comprehensive theoretical modelling and ob¬ 
servational constraints of DIG are crucial for extracting more in¬ 
formation from deeper IFU data. 


4.3 Ionisation parameter 

The ionisation parameter is quantified as the ionising photon flux 
through a unit area divided by the local number density of hydro¬ 
gen atoms. The ionisation parameter can be measured by taking 
the ratio of high ionisation to low ionisation species of the same 
atom. Using the strong lines available in this study, we measure 
the ionisation parameter using [Olll] A5007/[0ll] AA3726,3729 
(hereafter, [OlIl]/[Oll]). A theoretical calibration was first pre¬ 
sented in KD02 and a more user-friendly parametrisation is given 
in Kobulnicky & Kewley (2004, hereafter KK04). The latter is 
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Figure 5. Left and middle panels: Metallicity gradients of individual CALIFA galaxies measured using the two different abundance diagnostics (see Sec¬ 
tion 4.1). The straight lines indicate the best fits, and the dashed lines indicate ztlrr errors. The errors of the intercepts and the slopes are estimated from 
bootstrapping (see Section 4.4). Right panels: ionisation parameter as a function of radius. The ionisation parameters are derived using the [O III]/[0 II] diag¬ 
nostic (KK04; see Section 4.3). Each point in these plots corresponds to one IFU spaxel with significant (> Scr) detections on all the emission lines associated 
with the diagnostics. Spaxels contaminated by non-thermal excitation or dominated by DIG emission are rejected (see Section 4.2). The vertical dashed lines 
correspond to the radial cutoff within which the data are not considered in constraining the disk metallicity gradients. 


deproj. radius [kpcj 
5 10 15 


20 



deproj. radius[R25] 
^ deproj. radius [kpc^ 


20 



deproj. radiusIRjsl 


deproj. radius [kpc] 

5 10 15 


20 



deproj. radius [R25] 

deproj. radius [kpcl 
5 10 15 


20 



deproj. radius [R25j 


deproj. radius [kpc] 



0.0 0.2 0.4 0.6 0.8 1.0 
deproj. radius [R25 


1.2 1.4 


0 

8.2 r 

-- 8.0 ■ 
w 

E 7.8: 

S 7-6 ■ 

S 7.4 ■ 

f 7.2; 

7.0 ^ 


^ deproj. radius [kpc^ 


20 


0.0 0.2 0.4 0.6 0.8 1.0 1.2 
deproj. radius [R25 


1.4 


Figure 6. Same as Figure 5, but for the S12 galaxies. Each point corresponds to one Fill region extracted from the IFU data (see Section 2.3 and S12 for 
details). 


adopted in this paper. As emphasised in KD02, in addition to ion¬ 
isation parameter, the [O IIl]/[0 II] ratio also strongly depends on 
metallicity that needs to be known a priori. We adopt the metallic¬ 
ity from the N202 calibration to derive the ionisation parameter for 
each spectrum. 

4.4 Measuring metallicity gradients 

To derive metallicity gradients, we first convert the flux ratios to 
metallicities. For bulge-dominated galaxies or those with obvious 
bar structures in the CALIFA sample, we do not include data within 
(0.1-0.2) X i? 25 . At the centres of these generally high-mass sys¬ 
tems, the optical spectrum is dominated by the stellar continuum 
with strong B aimer absorption originating from an old stellar popu¬ 
lation. Correcting for Balmer absorptions in these regions is less ro¬ 


bust and more sensitive to both the SSP models and the algorithms 
adopted (e.g., Cid Fernandes et al. 2014, and references therein). 
The WiFeS galaxies do not suffer from this contamination because 
these low-mass systems are typically emission dominated even at 
the centres. The S12 and RIO galaxies also do not suffer from this 
artefact because the emission lines are directly measured towards 
Hll regions. 

To estimate the errors in the metallicity gradients, we 
adopt a bootstrapping approach similar to Kewley et al. (2010), 
Rupkeetal. (2010b) and Rich et al. (2012). We randomly draw 
from the measurements the same number of data points but with re¬ 
placement, and perform an unweighted least-squares linear fit with 
the drawn data. In each galaxy, this process is repeated 1,000 times 
and each fit result is recorded. The median and standard deviation 
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of the slopes and intercepts are considered as the best estimates of 
the metallicity gradient. 


5 RESULT 

5.1 Metallicity gradients of individual galaxies 

In Figure 5 and Figure 6, we present 4 metallicity gradients mea¬ 
sured in two CALIFA and two S12 galaxies, respectively. The two 
different metallicites are shown in the first two panels, and the last 
panels show radial profiles of the ionisation parameter. In the first 
two panels, straight lines indicate the best fits of the gradients, and 
dashed lines indicate Icr errors as propagated using analytic ex¬ 
pressions with bootstrapped errors. The rest of the CALIFA and 
S12 galaxies are presented in Figures Al and A2 in Appendix A. 

The metallicity gradients and the radial profiles of the ioni¬ 
sation parameter of our 10 WiFeS galaxies are also presented in 
Figure A3. The metallicity of the majority of WiFeS galaxies have 
metallicities 12-|-log(0/H) < 8.4 (in KD02 scale), where the con¬ 
version from 03N2 to N202 is difficult to determine, and therefore 
we are only able to derive their N202 metallicity gradients. These 
galaxies are excluded from the rest of the comparisons between two 
metallicity calibrations, but we include them later while comparing 
metallicity gradients of different galaxies (Section 5.4). 

In Figure 7, we compare metallicity gradients derived 
from 03N2 and from N202 for the CALIFA and S12 galax¬ 
ies. All the galaxies exhibit negative metallicity gradients in 
both calibrations, and 33%/73% of the galaxies agrees within 
±0.05/±0.14 dex These values can be considered as the 
level of residual systematics in the metallicity calibrations. Notice¬ 
ably there are several outliers well below the one-to-one line which 
we label in Figure 7. These objects could provide insight into the 
cause of the disagreement between the two calibrations. We fur¬ 
ther investigate the cause of this discrepancy in the following two 
subsections. 



N202 metallicity gradient [dex R 25 ] 


Figure 7. Left: Comparison between the metallicity gradients (dex ) 
derived using the 03N2 and the N202 diagnostics (see Section 4.1). The 
metallicity gradients of the majority of the galaxies (33%/73%) agree within 
±0.05/ib0.14 dex R^^ ■ The metallicity gradients of the four outliers la¬ 
beled are shown in Figure 5 and Figure 6. 


5.2 The effect of ionisation parameter 

In the left panel of Figure 8, we compare the metallicities derived 
using the 03N2 and N202 calibrations. Each point represents one 
measurement, i.e., one spaxel from the CALIFA galaxies or one Hll 
region from the S12 galaxies. We colour-code the points by logfg) 
derived using the [O lll]/[011] (KK04) calibration. The left panel 
of Figure 8 clearly demonstrates that superficially the two metal¬ 
licities agree reasonably well; 57%/83% of the metallicities agree 
within ±0.05/0.1 dex. Similar comparisons were also carried out 
by Rupke et al. (2010b) where comparable scattering also exists. 

In the left panel of Figure 8, the scattering around the one- 
to-one line is not random, but correlates with ionisation param¬ 
eter. Spaxels or slits with high ionisation parameter have higher 
N202 metallicities than 03N2 metallicities. Those with low ion¬ 
isation parameter have lower N202 metallicity than 03N2 metal¬ 
licities. In the right panel of Figure 8, we show the difference be¬ 
tween metallicities derived with the 03N2 and N202 ratios ver¬ 
sus ionisation parameter. The differences between the two diag¬ 
nostics can be up to 0.2 - 0.4 dex at extreme values of ionisa¬ 
tion parameter (log(g) < 7.0 cm or log(g) > 8.2 cm 
). At log]^) > (<)7.3 cm s“^, the 03N2 diagnostic gives higher 
(lower) metallicity than the N202 diagnostic. Clearly, ionisation 
parameter is the cause of the discrepancy. As we emphasised ear¬ 
lier, 03N2 is calibrated empirically without taken into account the 
change of ionisation parameter. A theoretical calibration of 03N2 
will be presented in Kewley et al. (in preparation) and will rec¬ 
oncile this discrepancy with new stellar population synthesis and 
photoionisation models. 


5.3 Discrepancies among metallicity gradients 

We now return to discuss the cause of the discrepancies in the 
metallicity gradients shown in Figure 7. 

Figure 5 (and Figure 6) show metallicity gradients for two 
CALIFA (two S12 galaxies) that are labeled as outliers in Figure 7. 
These galaxies have steeper 03N2 than N202 metallicity gradients 
because, at large radii, 03N2 metallicities are systematically lower 
than N202 metallicities. As shown in Figure 8, lower 03N2 than 
N202 metallicities naturally arises when log(g) > 7.3 cm s“^. 
Indeed in the third panels of Figure 5 and Figure 6, these galaxies 
typically have log(g) > 7.3 cm at large radii. These galax¬ 
ies all show indications of smooth rising of ionisation parameter 
from their centres to outskirts, implying a continuous radial change 
of their properties of the ionising radiation. The higher ionisation 
parameters could be caused by the more active star formation ac¬ 
tivities with different distributions of molecular gas (Dopita et al. 
2014). 

In extreme cases, the differences in metallicity gradient mea¬ 
sured with N202 and 03N2 can be up to a ~ 0.4 dex 
(e.g., NGC3184 in Figure 6). Similar findings are also reported 
in Rupke et al. (2010b). These results have important implica¬ 
tions for metallicity gradient studies at high redshift, where typ¬ 
ically only [N II] AA6548,6583 and Ha are available (in some 
cases also [Olll] AA4959,5007 and H/3, e.g., Crescietal. 2010; 
Yuan et al. 2011; Jones et al. 2010, 2013; Queyrel et al. 2012; 
Swinbank et al. 2012). While all the diagnostics using these four 
lines, i.e. [N II] A6583/Ha and 03N2, are sensitive to the change of 
ionisation parameter, quantitative interpretation of metallicity gra¬ 
dients should bear in mind the potential impact of ionisation pa¬ 
rameter gradients in galaxies. 

For galaxies not labeled in Figure 7, we do not find obvious 
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Figure 8. Left: Comparison between the metallicities derived using the 03N2 and N202 diagnostics (see Section 4.1). Each CALIFA data point is an IFU 
spaxel, and each S12 data point is an Hit region extracted from IFU data. Data points are colour-coded with their corresponding ionisation parameters (see 
Section 4.3). A total of 57%/83% of the data points agrees within ±0.05/0.1 dex. The degree of disagreement correlates strongly with the ionisation parameter. 
Right: Difference between the metallicities derived using the 03N2 and N202 diagnostics versus the ionisation parameter. A strong anti-con'elation between 
the two quantities is obvious. More discussion about this discrepancy of metallicities is provided in Section 5.2, and the impact on measuring metallicity 
gradients in Section 5.3. 


signs of a correlation between ionisation parameter and radius. The 
inconsistency in metallicities derived using the 03N2 and N202 
diagnostics does not correlate with radius. When measuring metal¬ 
licity gradients in these galaxies, the net effect is merely to increase 
the uncertainties in the measurements, rather than biasing the mea¬ 
surements in any systematic way. This is essentially the cause of 
the small scatter (< 0.1 dex in Figure 7. 

5.4 Metallicity gradients in field star-forming galaxies 

After investigating the systematics induced by the variation of ion¬ 
isation parameter, we adopt the N202 measurements as our final 
metallicity gradients and we now compare metallicity gradients in 
field star-forming galaxies with their stellar mass and 6-band lumi¬ 
nosity. 


5.4.1 Metallicity gradient - stellar mass 

The upper and lower panels of Figure 9 show the metallicity gradi¬ 
ent versus stellar mass of our four samples in units of dex kpc~^ 
and dex Riff respectively. In the upper panel, the metallicity gradi¬ 
ents appear to depend on stellar mass where (1) low mass galaxies 
have on average a steeper metallicity gradient, and (2) are more 
diverse in the steepness of metallicity gradients compared to high- 
mass galaxies. To quantify the dependency, we split the sample into 
two mass bins with roughly equal numbers of galaxies, i.e. a high- 
mass bin (log{Mt/M q) > 9.6; Ngai = 24) and a low-mass bin 
{log{Mt/M q) < 9.6; Ngai ~ 25). We compute the bootstrapped 
means and standard deviations of the metallicity gradients. The re¬ 
sults are tabulated in Table 1. We indeed find that the low-mass 
galaxies have a steeper mean metallicity gradient at 3.4a level than 
the high-mass galaxies, and also a larger standard deviation of the 
metallicity gradients at 3.7a level. 

When the metallicity gradients are normalised to the galaxy 
sizes (i.e. dex Rff, the lower panel of Figure 9 does not support 
any clear dependency of metallicity gradient on stellar mass. The 
difference between the means is only at 1.2cr level, with low-mass 
galaxies having slightly flatter metallicity gradients than high-mass 
galaxies. The standard deviations of the metallicity gradients are 
virtually identical (the difference is only 0.6 ct). 
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Figure 9. Metallicity gradient versus stellar mass when the metallicity gra¬ 
dients are measured in dex kpc“ ^ (upper panel) and in dex Rf (lower 
panel). More details are discussed in Section 5.4.1. 


The different dependency of metallicity gradient on stellar 
mass while measuring the metallicity gradient in absolute scale 
(kpc) or relative scale (R 25 ) can be understood as a size effect. If 
galaxies with steeper dex kpc“^ metallicity gradients are smaller 
in their physical sizes (small R 25 ), then the steep dex kpc“^ metal¬ 
licity gradients would be compensated when the galaxy sizes are 
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Figure 10. Metallicity gradients (dex kpc“^) versus R 25 in kpc. A corre- 
lation between the two quantities can be seen. More details are discussed in 
Section 5.4.1. 
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Figure 11. 6-band luminosity versus R 25 in kpc. The dashed line indicates 
the luminosity-size relation by van den Bergh (2008). 


taken into account. Figure 10 shows the dex kpc“^ metallicity 
gradient versus galaxy size R 25 of our samples. Indeed, galaxies 
with steeper metallicity gradients generally have smaller R 25 than 
galaxies with shallower metallicity gradients. The Spearman rank 
correlation coefficient between R 25 and the dex kpc“^ metallicity 
gradients is 0.6, which is different from zero (i.e. no correlation) at 
a significance of 7.5 x 10~®. In Figure 11 , we demonstrate that our 
samples fall on the luminosity-size relation (van den Bergh 2008), 
i.e. low luminosity galaxies have smaller R 25 than high luminosity 
galaxies, indicating that the steep dex kpc“^ metallicity gradients 
in low-mass galaxies could he associated with their small physical 
sizes. These results imply that the evolution of metallicity gradients 
is closely related to the growth of galaxy size. 


5.4.2 Metallicity gradient - Mb 

We now compare metallicity gradients with the absolute B- 
band magnitudes, Mb, for our four samples. The absolute B- 
band magnitude is used as a proxy for mass in metallicity stud- 
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Figure 12. Metallicity gradient versus absolute B-band magnitude when 
the metallicity gradients are measured in dex kpc“^ (upper panel) and in 
dex R^^ (lower panel). More details are discussed in Section 5.4.2. 


ies where multi-band photometry is not available. Rubin et al. 
(1984) first showed that metallicity is correlated with luminosity 
in disk galaxies. Further investigations solidified the luminosity- 
metallicity correlation in nearby disk galaxies (Bothun et al. 1984; 
Wyse & Silk 1985; Skillman et al. 1989; Vila-Costas & Edmunds 
1992; Zaritsky et al. 1994; Garnett 2002). We emphasise that opti¬ 
cal luminosity is not always a reliable surrogate for the stellar mass 
of a galaxy because optical luminosities are sensitive to the cur¬ 
rent SFR and are affected by dust. Near-infrared luminosities can 
be influenced by the age of the stellar population of a galaxy. Since 
this quantity is widely adopted in earlier studies, we present our 
measurements below for comparison. 

In Figure 12, we present metallicity gradient versus Mb of our 
four samples. Metallicity gradients are shown in both dex kpc“^ 
(upper panel) and dex R^^ (lower panel). Similarly, we split the 
sample into a high luminosity bin (Mb < —20.1; Ngai = 24) 
and a low luminosity bin (Mb > —20.1; Ngai = 25), and 
compute the bootstrapped means and standard deviations of the 
metallicity gradients (Table 1). We reach the similar conclusions, 
as in the comparisons with stellar mass, that (1) low luminos¬ 
ity galaxies have a steeper mean dex kpc” metallicity gradi¬ 
ent (3.4 (t), and (2) low luminosity galaxies have a larger stan¬ 
dard deviation of dex kpc”^ metallicity gradients (3.5cr). When 
the galaxy sizes are taken into account, i.e. dex R^^, the low and 
high luminosity galaxies have very similar mean metallicity gra¬ 
dients (1.4 ct) and standard deviations (1.8cr). These findings are 
in agreement with previous studies (Vila-Costas & Edmunds 1992; 
Garnett et al. 1997 ; Prantzos & Boissier 2000). The lack of correla- 
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Table 1. A local benchmark gradient 



Mean 

Standai'd deviation 

dex kpc ^ 

log(M*/M0) > 9.6 

-0.026 ± 0.002 

0.010 ± 0.001 

\og(Mt/M q) < 9.6 

-0.064 ±0.011 

0.054 ± 0.012 

Mb < -20.1 

-0.025 ± 0.002 

0.008 ± 0.001 

Mb > -20.1 

-0.063 ±0.011 

0.053 ±0.013 

dex -^ 25 ^ 

\og(Mt/M q) > 9.6 

-0.42 ±0.03 

0.16 ±0.03 

\og{Mt/M q) < 9.6 

-0.36 ±0.04 

0.18 ±0.02 

Mb < -20.1 

-0.40 ±0.03 

0.12 ±0.02 

Mb > -20.1 

-0.34 ±0.03 

0.17 ±0.02 

AlH 

-0.39“ 

0.18“ 


Means, standard deviations and the associated en'ors are derived from boot¬ 
strapping. 

‘‘Mean and standard deviation of all the galaxies (i.e. not from bootstrap¬ 
ping). See Figure 13. 

lions between the metallicity gradient (in dex per scale-length) and 
macroscopic properties such as stellar mass and 6-band magnitude 
might imply that the relationships between these parameters are 
more complex than simple correlations, a hypothesis first suggested 
by Zaritsky et al. (1994) and recently investigated by Pilyugin et al. 
(2014b). 

Few et al. (2012) performed cosmological zoom-in simula¬ 
tions that yield metallicity gradients, disk scale-lengths, 6-band 
magnitudes, and stellar masses of 19 galaxies in field and loose 
group environments. The simulations focused on Milky Way-mass 
galaxies spanning stellar mass and 6-band magnitude ranges of 
10.4 < log(M*/Mo) < 11.1 and -19.7 > Mb > -21.7, re¬ 
spectively, which correspond to the high-mass and high-luminosity 
ends of our samples. Few et al. (2012) found no significant dif¬ 
ferences in metallicity gradients between the galaxies in field 
and loose group. The overall metallicity gradients are —0.046 ± 
0.013 dex kpc“^ and —0.40 ± 0.13 dex (mean ± standard 
deviation). Here, we convert the disk scale lengths Rd reported 
by Few et al. (2012) to 7^25 assuming that the exponential disks 
have the canonical central surface brightness for normal spirals of 
21.65 magnitude/arcsec^ (Freeman 1970). The overall metallicity 
gradients by Few et al. (2012) are consistent with our results in the 
high-mass and high-luminosity ends. Similar simulations in the fu¬ 
ture targeting lower masses and luminosities could provide valuable 
constraints for the different prescriptions built into the simulations. 


6 A LOCAL BENCHMARK GRADIENT 

We provide a local benchmark of metallicity gradients inspired by 
the uniformity of the dex R^^ metallicity gradients. We believe 
the local benchmark gradient will be useful for comparison with 
metallicity gradients measured at high redshift. In Figure 13, we 
present the distribution of the 49 measured metallicity gradients; 
and we summarise the mean and standard deviation in Table 1 . The 
benchmark metallicity gradient measures —0.39 ± 0.18 dex R^^ 
(mean ± standard deviation). A one-sided Kolmogorov-Smimov 
test yields a 78% probability for the observed distribution to be 
drawn from the normal distribution shown as the black curve in Fig¬ 
ure 13. We note that the difference between the 03N2 and N202 
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Figure 13. Distribution of the 49 metallicity gradients. The overall 
mean and standard deviation of the metallicity gradients are —0.39 ± 
0.18 dex 6 ^ 5 ^ (i.e. the benchmark metallicity gradient). The black curve 
indicates a Gaussian with these characteristic values, i.e. not a fit to the 
distribution. A one-sided Kolmogorov-Smirnov test yields a probability of 
78% for the observed distribution to be drawn from the back curve. 


metallicity gradients is typically 0.14 dex R^^ (Figure 7), compa¬ 
rable to the standard deviation of 0.18 dex R^^ in the benchmark 
gradient. This remarkably small difference suggests that the intrin¬ 
sic spread of the metallicity gradients could be even tighter than 
0.18 dex R 2 ^ although precisely quantifying the tightness is non¬ 
trivial due to the various systematic effects. 

We note that in Figures 9 and 12, the RIO sample appears to 
have lower metallicity gradients than the other three samples. A 
one-sided Kolmogorov-Smimov test comparing the RIO sample to 
the benchmark gradient, and a two-sided Kolmogorov-Smirnov test 
comparing the RIO sample to the whole sample both reveal moder¬ 
ately low but non-zero probability (10%) for the two distributions 
to be the same. It is possible that the discrepancy is due to low 
number statistics. Furthermore, a two-sided Kolmogorov-Smimov 
test comparing the full sample to that with the RIO galaxies ex¬ 
cluded yields a high probability (99%) for the two distributions to 
be drawn from the same parent distribution. We conclude that the 
discrepancy is only apparent and does not change our results. 

S12 measured metallicity gradients in 25 face-on spirals and 
found a common metallicity gradient of —0.12 ±0.11 dex RJ^ 
(median ± standard deviation). Sanchez et al. (2014) expanded 
the study to 193 galaxies with the CALIFA survey, and found a very 
similar common metallicity gradient of —0.10 ± 0.09 dex R^^- 
Sanchez et al. (2014) showed that, with their large sample size, 
the metallicity gradients are independent of morphology, incidence 
of bars, absolute magnitude and mass, a result that had also been 
hinted by S12. Sanchez et al. (2014) found that the only clear cor¬ 
relation is between merger stage and metallicity gradient, where the 
slope is flattened as merger progresses (see also Kewley et al. 2010; 
Rupke et al. 2008, 2010a,b; Torrey et al. 2012; Rich et al. 2012). 

To compare the common metallicity gradient by Sanchez et al. 
(2014) with our benchmark gradient, one must take into ac¬ 
count the different scale-lengths, i.e. R^ and R 25 , and metallic¬ 
ity calibrations (Sanchez et al. 2014 and S12 both adopted the 
O3N2/PP04 calibration). Assuming again an exponential disk with 
the canonical central surface brightness for normal spirals (i.e. 
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Figure 14. Comparison between our benchmark gradient (dashed curve) 
and the metallicity gradients from 104 field spiral galaxies published by 
Pilyugin et al. (2014a). Our benchmark gradient is —0.39 zb 0.18 dex 
(mean ± standard deviation; Figure 13). The metallicity gradients of 
the field spiral galaxies from Pilyugin et al. (2014a) measure —0.32 di 
0.20 dexfl-l. 


21.65 magnitude/arcsec^; Freeman 1970) and the empirical con¬ 
version between the O3N2/PP04 and N2O2/KD02 metallicities 
(Kewley & Ellison 2008), we can convert the common metallicity 
gradient by Sanchez et al. (2014) to —0.20 ±0.18 dex , which 
is about a factor of two shallower than our benchmark gradient. 
The difference could be caused by the presence of close pairs and 
mergers in their sample. In addition, systematic errors such as the 
metallicity calibrations, variations of the ionisation parameter, and 
flattening of metallicity gradient at larger radii could all affect the 
measured metallicity gradients. 

Pilyugin et al. (2014a) complied more than 3,000 published 
spectra of Hit regions and adopted the calibration proposed by 
Pilyugin et al. (2012) to measure the metallicity gradients of 130 
nearby late-type galaxies. Using 104 of their field spiral galax¬ 
ies (i.e. excluding mergers and close pairs), we constrain a mean 
and standard deviation of the metallicity gradients of —0.32 ± 
0.20 dex ■ These values are consistent with our benchmark 
gradient (see Figure 14 for a comparison), but the distribution from 
their sample do not match our benchmark gradient, implying resid¬ 
ual systematic errors. 


7 WHY IS THERE A BENCHMARK GRADIENT? 

The existence of a common metallicity gradient when expressed 
with respect to some scale-lengths had long been suggested (e.g., 
Zaritsky et al. 1994; Garnett etal. 1997; Vila-Costas & Edmunds 
1992), and is further solidified by recent observations with inte¬ 
gral field spectroscopy on large samples (S12; Sanchez et al. 2014). 
The common slope implies that all disk galaxies went through 
very similar chemical evolution when building up their disks, pre¬ 
sumably in an inside-out fashion (Sanchez et al. 2014). We show 
that our benchmark gradient is closely related to the growth of 
galaxy size, indicating that the chemical richness of disk galaxies 
co-evolves with the increase in their spatial dimensions (ses also 
Prantzos & Boissier 2000). 

Since the measured metallicity is the ratio of oxygen to hydro¬ 
gen atoms, to the zeroth-order the metallicity ought to be a strong 


function of the stellar-to-gas mass ratio. The stellar mass traces the 
total amount of metals produced through the formation of stars that 
drive nucleosynthesis; and the gas mass serves as the normalisation 
for the definition of metallicity. In the simplest case, known as the 
“closed-box” model, the chemical evolution began with pristine gas 
and experienced no subsequent mass exchange with the material 
outside the box. This classical closed-box model (Searle & Sargent 
1972; Pagel & Patchett 1975) already encapsulates the close link 
between the metallicity, Z, and the observed stellar mass to gas 
mass ratio, M^o/Mg\ 


m 


y 

l-R 


In 


. Mgit) 


y 

l-R 


In 


1 


(5) 

(6) 


Here, Z is the mass ratio instead of the number ratio adopted in 
12 ± log(0/H), y is the nucleosynthesis yield and R is the stellar 
returned mass fraction. The “observed” stellar mass, M*o, takes 
into account the mass loss through stellar winds described by the 
returned mass fraction, i.e. M*o = (1 — R)M*, where the time 
derivative of M* is the star formation rate. The gas fraction y,g is 
defined as 


a (t) = 

Such simple picture, however, is usually complicated by gas in¬ 
flows (e.g., accretion and merger) and outflows (e.g., AGN and star- 
burst feedback) that can remove or replenish both gas and metals to 
the ISM and break down the closed-box assumption. 

It is possible to put constraints on the inflow and out¬ 
flow history by studying the global metallicity of star-forming 
galaxies using the mass-metallicity relation (e.g., Spitoni et al. 
2010; Peeples & Shankar 2011; Zahid et al. 2012; Lilly et al. 2013; 
Peeples et al. 2014; see also Dave et al. 2011a). Global studies of 
metallicity suggest a close relation between metallicity, gas mass, 
and stellar mass, and the potential of the co-evolution of these three 
quantities following a simple, universal relation between metallic¬ 
ity and stellar-to-gas mass ratio (Zahid et al. 2014c; Ascasibar et al. 
2014). Similar universal relations may also exist on spatially re¬ 
solved scales (Ascasibar et al. 2014; see also Rosales-Ortega et al. 
2012 ). 

If such relation exists, the common metallicity gradient im¬ 
plies a close link between the radial profiles of the stars and the 
gas. Indeed, carbon monoxide and Hi 21 cm observations in nearby 
spiral galaxies reveal that the neutral (molecular and atomic) gas 
surface density profiles, Sg(»'), exhibit a tight universal profile 
(Bigiel & Blitz 2012). When the gas surface density profiles of the 
individual galaxies are expressed in terms of R 25 and normalised 
to a transition radius where the molecular and atomic gas have the 
same surface densities, Bigiel & Blitz (2012) show that the overall 
surface density of all the galaxies (at r > 0 . 2 i? 25 ) follows a sim¬ 
ple exponential profile. The exponential profile has a logarithmic 
slope of —0.71 dex R^^ and a very small bootstrapped error of the 
mean of 0.06 dex R^^ (see their figure 3). In Figure 15, we show 
the gas, the stellar (E*(r)), and the stellar-to-gas surface density 
profiles of 14 field spiral galaxies from Leroy et al. (2008, 9.3 < 
log(M,/M q) < 10.9; R 25 = 3-20 kpc; c.f. Figures 9 and 10). 
All the profiles are expressed in terms of R 25 , and we normalise the 
gas and stellar surface density profiles at 0.6i?25. Clearly, the gas, 
stellar, and stellar-to-gas surface density all follow simple radial 
profiles with common slopes. We fit simple exponential profiles 
to each galaxy using data at r > 0 . 2 i ?25 to constrain the loga- 
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Figure 15. Gas (atomic + molecular; top panel), stellar (middle panel), and 
stellar-to-gas (bottom panel) surface density profiles of 14 field spiral galax¬ 
ies from Leroy et al. (2008). The radial distance is expressed in terms of 
U 25 . The gas and stellar profiles are normalised at O. 6 R 25 , and the stellar- 
to-gas profiles are not normalised. The dashed line in the top panel indicates 
the best-fit universal gas profile from Bigiel & Blitz (2012). The insets show 
the distributions of the slopes measured by fitting straight lines to the loga¬ 
rithmic profiles using data at r > 0 . 2 i? 25 . 


rithmic slopes, and we find that the distributions of the slopes have 
means ± standard deviations of —0.69 ±0.40 dex for the gas, 
— 1.68 ± 0.23 dex for the stars, and —0.98 ± 0.35 dex 
for the stellar-to-gas surface density. 

These tight universal profiles, in particular the stellar-to-gas 
mass ratio, may govern the common metallicity gradient in disk 
galaxies. We use simple chemical evolution models to quantita¬ 
tively address the relation between metallicity, and stellar-to-gas 
mass ratio. The models we adopted are special cases of the more 
general derivations by Recchi et al. (2008). Similar models have 
been applied to global metallicities of galaxies (e.g., Spitoni et al. 
2010; Dayal et al. 2013; Lilly et al. 2013; Pipino et al. 2014). The 
analytical models consider, for a given volume element, metal pro¬ 
duction by stars, stellar mass return, inflows, and outflows, under 
the assumption that the gas is well-mixed and the mass return from 
stars is recycled instantaneously. The inflows and outflows are de¬ 
scribed through two critical parameters: the mass loading factor 


n = 


Mlos. 

%jl 


( 8 ) 
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Figure 16. Examples of our chemical evolution models. The upper panel 
shows the r]-A plane. Each colour point con'esponds to one model deter¬ 
mined by the set of mass loading and mass accretion factors. The lower 
panel show the coiresponding relationship between the metallicity and 
stellar-to-gas mass ratio. The black dot at (rj, A) = (0, 0) in the upper 
panel is the closed-box model; the corresponding curve in the lower panel 
is shown as the dashed curve. The gray band marks the ±1 (t range of the 73 
massive star-forming galaxies measured by Tacconi et al. (2013). As only 
the molecular gas was measured, not the atomic gas, the range represents 
an upper limit. We note that denotes the “observed” stellar mass taken 
into account the stellar mass return. 


and the mass accretion factor 


^rtaccr 

ijj 


(9) 


where Miosa and Maccr are the mass-loss and mass-gain rates, re¬ 
spectively, and Ip is the SFR. Two additional assumptions are made. 
First, we assume that the accreted gas is metal-free and the outflow¬ 
ing gas has the same metallicity as the ISM at the time the outflows 
are launched. Second, we require the mass loading and mass ac¬ 
cretion factors to be constant. Note that this assumption only con¬ 
straints the factors (ratios) to be constant, and does not restrict the 
time evolution of star formation, inflow rate and outflow rate to any 
specific forms. 

The first assumption of outflows is valid if most of the out¬ 
flowing gas is entrained ISM close to the energy sources (e.g. super- 
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novae). Indeed, more than 75% of the outflowing gas is estimated 
to be entrained gas in nearby mergers (Rupke & Veilleux 2013). 
The second assumption of constant mass-loading and mass accre¬ 
tion factors helps the models to remain analytic and simple, and 
the assumption carries significant physical meanings. The constant 
mass-loading factor is postulated because the energy driving the 
outflowing mass is from star formation, and rj simply reflects the ef¬ 
ficiency in transferring energy from star formation to the outflowing 
gas. Similarly, because the inflowing gas supplies the reservoir for 
star formation, the constant mass accretion factor can be realised as 
the efficiency of collapsing the gas into stars. We explore different 
rj and A values later to understand the possible impact of them not 
being constant over time. Recchi et al. (2008) explore non-constant 
A by assuming exponential inflows and a linear Schmidt law. They 
conclude that constant A is a reasonable approximation of the late 
evolution of a galaxy provided that the infall timescale is of the 
same order of magnitude of the star formation timescale. 

Under these assumptions, a pair of non-negative (jj, A) deter¬ 
mines a unique analytical solution for the metallicity and the stellar- 
to-gas mass ratio: 


Z{t) = 



1-f 



t?-A \ 

1-r) 


where 


M,{t) 



V>0, 


A > 0, and 


rj — A 
1-R 


7 ^- 1 . 


( 11 ) 


That is, the stellar-to-gas mass ratio at any given time dictates the 
metallicity at that instant. The reader is referred to Kudritzki et al. 
(in preparation) for derivation of the models and special cases when 
the conditions in Equation 11 are not met. We also show in Ku¬ 
dritzki et al. (in preparation) that in the trivial case where there are 
no inflows nor outflows, i.e. rj — 0 and A = 0, our model is identi¬ 
cal to the classical closed-box model (i.e. Equation 5). 

In Figure 16, we present the relationships between the metal¬ 
licity and stellar-to-gas mass ratio (lower panel) using different sets 
of {rj, A) values (upper panel). We adopt the yield for the oxygen of 
0.00313 and the returned mass fraction of 0.4, and we will discuss 
the systematics of these two constants later. The classical closed- 
box model, which corresponds to the origin on the rj - A plane 
(upper panel), is also shown as the dashed line in the lower panel. 
In the gas rich regime, i.e. log(M,o/Mg) <C —0.5, all the mod¬ 
els coalesce as, under vast gas reservoirs, inflows and outflows do 
not change the metallicity appreciably. In the gas poor regime, i.e. 
\og{Mto/Mg) 3> —0.5, for a given stellar-to-gas mass ratio the 
metallicity is sensitive to the adopted mass loading and mass ac¬ 
cretion factors. This regime can be explored with the the 14 field 
spiral galaxies from Leroy et al. (2008) (0 < log(E*/Et,) < 2; c.f. 
Figure 15), which could place constraints on the mass loading and 
mass accretion factors. 

By combining the models and the measured stellar-to-gas 
mass profiles of the 14 field spiral galaxies, we can predict their 
metallicity gradients. The predicted metallicity gradients can be 
compared with our benchmark gradient to place constraints on the 
models. We calculate the predicted metallicity gradients by first 
converting the stellar-to-gas mass ratios to metallicities for each 
radial bins, and we fit linear profiles to data at r > 0.27^25 in each 
galaxies to derive the predicted metallicity gradients. In Figure 17, 
we show the distributions of the predicted metallicity gradients us¬ 
ing a grid of irj,A). In the top and middle panel of Figure 18, 
we compare the means and standard deviations of the predicted 


metallicity gradients to those from our benchmark gradient. In the 
bottom panel, we perform one-sample Kolmogorov-Smimov tests 
to compare the distributions of the predicted metallicity gradients 
with the benchmark gradient that has a Gaussian distribution of 
—0.39 ± 0.18 dex R^^ (mean ± standard deviation; Figure 13). 

Figure 17 qualitatively demonstrates that small rj and A val¬ 
ues (panels toward the lower left comer) are preferred because 
the predicted distributions are similar to the benchmark gradient. 
Large rj and A values tend to overproduce flatter metallicity gra¬ 
dients, effectively shifting the means of the distributions towards 
zero and reducing the widths of the distributions. This behaviour 
can be trivially understood with the bottom panel of Figure 16 
where the models with large rj and A values flatten at high stellar- 
to-gas mass ratios, yielding the same metallicities across the disks, 
i.e. flat metallicity gradients. We quantitatively address the allowed 
rj and A values in Figure 18 by investigating the differences of 
the means (upper panel), those of the standard deviations (middle 
panel), and the probabilities of reproducing the benchmark gradient 
through the Kolmogorov-Smirnov test (bottom panel). We find that 
for 0 < A < 0.2 and 0 < rj < 2, the Kolmogorov-Smimov tests 
yield good probabilities (> 20%) for the distributions of the pre¬ 
dicted metallicity gradients to be drawn from the benchmark gra¬ 
dient. Within the same ranges, the differences between the mean 
of the benchmark gradient and those of the predicted gradients are 
within about 20%, and the differences between the standard devi¬ 
ations are also within about 20%. Interestingly, the differences of 
the means, and those of the standard deviations both show that the 
closed-box model is the best model, but the Kolmogorov-Smimov 
tests suggest that low (but non-zero) mass loading and mass accre¬ 
tion factors are more preferred. While the precise values of rj and 
A probably cannot be determined from the 14 galaxies alone (due 
to low number statistics) and are likely to vary from system to sys¬ 
tem, closed-box and virtually closed-box are the models that can 
successfully reproduce the observed benchmark gradient. 

The success of reproducing the benchmark gradient with our 
simple models, however, do not imply that all galaxies evolve as 
closed-box or virtually closed-box throughout their lifetime. Nor 
do our results support the idea that galaxies always have con¬ 
stant mass loading and mass accretion factors, either on global 
or spatially resolved scales. Observations of high-redshift galax¬ 
ies (z > 1) have provided evidences that galaxies in the early 
Universe undergo many, perhaps intermittent, accretion events, 
immense star formation and outflows (e.g., Tacconi et al. 2010, 
2013; Weiner etal. 2009; Steidel et al. 2010; Genzel et al. 2011; 
Newman et al. 2012; Genzel et al. 2014). Similar outflows, in par¬ 
ticular the starburst-driven winds, are found to be ubiquitously in 
galaxies at lower redshifts with high enough star-formation sur¬ 
face density (Heckman 2002), with the wind velocities showing 
indication of correlating with both the SFR and host galaxy mass 
(Rupke et al. 2005a,b; Veilleux et al. 2005; Chen et al. 2010). En¬ 
ergy and mass return from both the stars and AGNs (i.e. “feed¬ 
back”) are indispensable for numerical simulations to reproduce 
many observed properties of galaxies, such as the stellar mass 
function and mass-metallicity relation (e.g., Springel & Hernquist 
2003; Oppenheimer et al. 2010; Dave et al. 201 la, b). From the the¬ 
oretical considerations, Murray et al. (2005) suggest that the mass 
loading factor could vary with the host galaxy mass, following 
different scaling relations depending on the winds being energy 
or momentum-driven. The later is favoured by recent smoothed- 
particle hydrodynamics -l- N-body simulations on both galactic 
and cosmological scales (e.g., Dave et al. 2011a,b; Hopkins et al. 
2012 ), but the two mechanisms dominating in galaxies of different 
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Figure 17. Distributions of the predicted metallicity gradients of the 14 field spiral galaxies shown in Figure 15. The metallicity gradients are predicted using 
the analytical models described in the text and in Kudritzki et al. (in preparation). A grid of mass loading factor rj and mass accretion factor A is adopted to 
predict the metallicity gradients. The r} and A values are labeled in the outer, large axes. The closed-box model is at (rj, A) = (0, 0). The measured benchmark 
gradient is shown as the black dashed curves for comparison. All the panels, except for the two top-left panels, have the same scales, as indicated in the 
bottom-left panel; the y-scales of the two top-left panels are labeled separately to accommodate the concentration of flat metallicity gradients in one bin. 


masses has also been suggested (Dutton & van den Bosch 2009). 
Unfortunately, measuring the mass loading factor accurately from 
observations remains difficult and a consensus on its values has not 
been reached yet (Zahid et al. 2014a). The multi-phase nature of 
the outflowing gas that spans wide ranges in both density and tem¬ 
perature poses a major observational challenge (see Veilleux et al. 
2005 for a review). 

Despite the complexity and the lack of observational con¬ 
straints on outflows and inflows, our simple models still can repro¬ 
duce the benchmark gradient because the metallicity is not sensi¬ 
tive to the adopted mass loading and accretion factors in the gas 
rich regime (log(M*o/Mg) <C —0.5; Figure 16). Recent radio 
observations of 2 > 1 galaxies reveal that galaxies at high red- 
shifts are typically gas rich (Tacconi et al. 2010, 2013), and there¬ 
fore the potentially high mass loading and mass accretion factors at 
high redshifts do not determine the metallicity gradients at 2 = 0. 
We hypothesise that as galaxies evolve to higher stellar-to-gas ratio 
(presumably at 2 < 1), both the mass loading and accretion factors 
decrease dramatically and stabilise such that their chemical evolu¬ 
tion can be approximated by the closed-box or virtually closed-box 
models. 

Our analysis favours a very low mass accretion factor (A < 
0.3), consistent with the lack of direct observational evidence of 
gas accretion in field galaxies in the local Universe. The low mass 
accretion factor and high stellar-to-gas mass ratio imply that field 


star-forming galaxies in the local Universe have no significant, re¬ 
cent refuelling of their gas reservoirs and their low level of star 
formation activities are sustained by the remaining gas reservoirs 
acquired presumably at high redshift. We also obtain a marginally 
low mass loading factor of about r] <2. Such mass loading factor 
is consistent with the range of rj measured. Zahid et al. (2012) em¬ 
pirically constrain the mass loading factor in star-forming galax¬ 
ies to be less than 1 by assuming that these galaxies evolve 
on the measured mass-metallicity and the galaxy main-sequence 
(Noeske et al. 2007). Bolatto et al. (2013) estimate the mass load¬ 
ing factor (of the molecular gas) of more than 1 (and probably 
~ 3) in the nearby starburst galaxies NGC 253. Mass loading 
factors of about 0.1-1 were also found in 2 < 0.5 luminous and 
ultra-luminous infrared galaxies, and nearby mergers (Rupke et al. 
2005b; Rupke & Veilleux 2013). 

We note that although we assumed the oxygen yield as con¬ 
stant, the oxygen yield varies with both the stellar metallicity 
and initial mass function (e.g., Maeder 1992; Woosley & Weaver 
1995; Kobayashi et al. 2006; see Zahid et al. 2012 for a summary). 
We also assumed a constant returned mass fraction, but the re¬ 
turned mass fraction is functions of both the stellar age and ini¬ 
tial mass function, spanning a range of approximately 0.15-0.45 
(e.g., Leitner & Kravtsov 2011), and 0.3-0.45 for the Salpeter and 
Chabrier initial mass function (Salpeter 1955; Chabrier 2003). Con¬ 
straining these two parameters individually has proven to be diffi- 
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Figure 18. A quantitative comparison between the benchmark gradient 
and the metallicity gradients predicted using the 14 field spiral galaxies 
from Leroy et al. (2008) and our analytical models. A qualitative compar¬ 
ison is also presented in 17. Each location on the plots corresponds to 
adopting one set of (rj, A) values to predict the 14 metallicity gradients. 
The top panel shows the differences in the means, and the middle panel 
shows the differences in the standard deviation, with contours indicating 
the differences in percentage. The bottom panel shows the p-values from 
one-sample Kolmogorov-Smirnov tests. A higher p-value indicates a higher 
probability for the distribution of the 14 predicted metallicity gradients to 
be drawn from the benchmark gradient that has a Gaussian distribution of 
—0.39 ± 0.18 dex (mean ± standard deviation). 


cult as they are degenerate through y/{l — R), i.e. the pre-factor 
in Equation 5 (see also Zahid et al. 2012). A higher degree of nu¬ 
cleosynthesis of the oxygen (higher y) can be balanced by lock¬ 
ing up more oxygen in each generation of stars (higher 1 — R; 
lower R), effectively leaving the same amount of oxygen in the 
ISM. In this work, we adopt the oxygen yield and return gas frac¬ 
tion from Kudritzki et al. (in preparation) who empirically constrain 
t//(l — i?) to the accuracy of 25% by reproducing the metallicity 
and the metallicity gradient of the young stellar population in the 
Milky Way. We vary the returned mass fraction from 0.15 to 0.45 
while keeping y/{l — R) fixed, which corresponds to an oxygen 
yield between 0.0045 to 0.0029, and our results do not change con¬ 
siderably. We find that a higher (lower) returned mass fraction re¬ 
sulting in more (less) gas return would flatten (steepen) the model 
curves in the bottom panel of Figure 16 at high stellar-to-gas mass 
ratio. However, the degree of flattening (steepening) is insignifi¬ 
cant such that similar preferred rj and A are recovered. For reason¬ 
able returned mass fractions {R = [0.45, 0.3, 0.15]), our preferred 


rj and A ranges, defined by the 20% contour of the Kolmogorov- 
Smirnov tests, remain virtually the same (0 < A < [0.2, 0.25, 0.3]; 

0<t7< [1.8,2.1,2.8]). 

In the models, we assumed that the outflowing gas has the 
same metallicity as the ISM at the time the outflows are launched. 
This assumption is appropriate because more than 75% of the out¬ 
flowing mass in nearby mergers is entrained gas (Rupke & Veilleux 
2013). Evidences of the hot, wind fluid being more enriched than 
the ISM have been reported at least in one nearby dwarf galaxy 
NGC 1569 (Martin et al. 2002), perhaps indicating that the hot ma¬ 
terials can survive the gravitational potential better than the cold 
entrained gas. Constraints on the metallicity of the outflowing gas 
remain scarce because X-ray observations are often required. If in¬ 
deed the outflowing gas has a higher metallicity than the ISM, the 
metallicity at a given stellar-to-mass ratio would be lower than that 
without a higher metallicity, particularly when the stellar-to-mass 
ratio is high. More enriched galactic winds, similarly, will flatten 
the model curves in the bottom panel of Eigure 16, causing our 
analysis to favour an even smaller mass loading factor. 


8 SUMMARY AND CONCLUSIONS 

We have presented metallicity gradients of 49 local field star¬ 
forming galaxies measured with integral field spectroscopy and slit 
spectroscopy. Metallicities have been determined for these galax¬ 
ies using strong optical emission lines ([Oil] AA3726,3729, H/S, 
[Out] A5007, Ha, and [Nil] A6583) with two widely adopted 
metallicity calibrations (the 03N2 diagnostic by Pettini & Pagel 
2004; and the N202 diagnostic by Kewley & Dopita 2002). Our re¬ 
sults show that the metallicities measured with the two calibrations 
are typically in good agreement (±0.1 dex), but the differences 
in metallicities correlate with the ionisation parameters. Similarly, 
the two calibrations yields metallicity gradients typically in good 
agreement (±0.14 dex R^^), but up to 0.4 dex R^^ difference 
is possible when the ionisation parameters change systematically 
with radius. 

When comparing the metallicity gradients with the stellar 
masses and absolute B-band magnitudes, we find that, when the 
metallicity gradients are expressed in dex kpc“^, galaxies with 
lower masses and luminosities have (1) on average a steeper metal¬ 
licity gradient and (2) more diverse metallicity gradients compared 
to galaxies of higher masses and luminosities. Such dependencies 
on mass and luminosity do not exist when the sizes of galaxies 
are taken into account and the metallicity gradients are expressed 
in terms of dex R^^. All our disk galaxies appear to have a com¬ 
mon metallicity gradient when normalised to the optical radii of the 
galaxies, consistent with previous studies. This leads us to quantify 
a local benchmark gradient of —0.39 ±0.18 dex R^^ that could be 
useful for comparison with metallicity gradients measured at high 
redshifts. 

We adopt simple chemical evolution models to investigate the 
cause of the common, uniform metallicity gradients. Starting from 
the measured atomic and molecular gas, and stellar surface density 
profiles in 14 nearby, field spiral galaxies, our analytical models 
can qualitatively and quantitatively reproduce the measured local 
benchmark gradient. Our results suggest that the galactic disks of 
spiral galaxies (at 0.2 < r/i ?25 < 1) evolve chemically close to 
the closed-box model when the stellar-to-gas ratio becomes high 
(log(M*o/Afg) S> —0.5). The inferred negligible mass accretion 
rates (< 0.3 x SFR), and very low mass outflow rates (< 3 x SFR) 
are broadly consistent with observational constraints. 
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To summarise, our simple chemical models already capture 
the fundamental physics governing the common metallicity gradi¬ 
ent. The common metallicity gradient is a direct result of the com¬ 
mon gas and stellar surface density profiles under the coevolution 
of gas, stars, and metals during galaxies build up their mass. 
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APPENDIX A: METALLICITY GRADIENTS AND 
IONISATION PARAMETERS OE INDIVIDUAL GALAXIES 

Figures A1 and A2 show the metallicity gradients and ionisation 
parameters for the rest of the CALIFA and S12 samples, respec¬ 
tively. These are in addition to the four galaxies presented in Fig¬ 
ures 5 and 6. Figure A3 show those for the WiFeS galaxies. 
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Figure Al. Continuation of Figure 5. 
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Figure A2. Continuation of Figure 6. 
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Figure A3. Same as Figure 5, but for the WiFeS galaxies. Only the N202 metallicity gradients (left panels) and ionisation parameter versus radius (right 
panels) are shown. See more details in Section 4.1. 
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